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ABSTRACT 


This  proposed  study  examines  the  potential  use  of  satellite  passive  microwave 
rainfall  measurements  derived  from  Special  Sensor  Microwave/Imager  (SSM/I) 
radiometers  onboard  the  Defense  Meteorological  Satellite  Program  (DMSP)  constellation 
to  improve  eastern  North  Pacific  Ocean  tropical  cyclone  intensity  change  forecasting 
techniques.  Relationships  between  parameters  obtained  from  an  operational  SSM/I-based 
rainfall  measuring  algorithm  and  12-,  24-,  36-,  48-,  60-  and  72-hour  intensity  changes 
from  best  track  data  records  are  examined  in  an  effort  to  identify  statistically  significant 
predictors  of  intensity  change. 

Correlations  between  rainfall  parameters  and  intensity  change  are  analyzed  using 
tropical  cyclone  data  from  three  years,  1992  to  1994.  Stratifications  based  upon  tropical 
cyclone  intensity,  rate  of  intensity  change,  climatology,  translation,  landfall  and  synoptic- 
scale  environmental  forcing  variables  are  studied  to  understand  factors  that  may  affect  a 
statistical  relationship  between  rainfall  parameters  and  intensity  change.  The  predictive 
skill  of  statistically  significant  rainfall  parameters  is  assessed  by  using  independent 
tropical  cyclone  data  from  another  year,  1995.  In  addition,  case  studies  on  individual 
tropical  cyclones  are  conducted  to  gain  insight  on  predictive  performance  and  operational 
implementation  issues. 
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Impacts  of  these  statistically  significant  rainfall  parameters  on  a  multiple  linear 
regression-based,  eastern  North  Pacific  tropical  cyclone  intensity  change  forecasting 
method  under  development  by  Hobgood  and  Petty  at  The  Ohio  State  University  are 
studied.  The  overall  goal  is  to  determine  if  SSM/I  rainfall  parameters  can  add  predictive 
skill  to  an  objective,  intensity  change  guidance  product.  The  guidance  product's  skill  in 
predicting  intensity  change  is  assessed  by  statistically  comparing  its  12-,  24-,  36-,  48-, 
60-  and  72-hour  intensity  change  forecast  errors  with  those  from  official  National 
Hurricane  Center  forecasts,  simple  persistence  and  no  change. 
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CHAPTER  1 


INTRODUCTION 


1.1  Defining  the  Problem 

Predictions  of  tropical  cyclone  movement  and  intensity  are  two  “critical”  forecast 
problems  (Simpson  and  Riehl  1981).  Currently,  intensity  forecasting  is  the  operational 
tropical  cyclone  forecaster’s  “greatest  problem”  requiring  the  development  of  effective 
objective  guidance  (Gross  and  Lawrence  1996).  This  research  proposes  to  examine  this 
important  problem  with  an  emphasis  on  improving  tropical  cyclone  intensity  prediction 
in  the  eastern  North  Pacific  Ocean  basin. 

Before  going  further,  some  of  the  relevant  terms  should  be  defined.  A  tropical 
cyclone  is  a  warm-core,  nonfrontal  low  pressure  system  of  synoptic  scale  that  develops 
over  tropical  or  subtropical  waters  and  has  a  definite  organized  surface  circulation 
(OFCM  1996).  This  study’s  basin  of  interest  is  the  Pacific  Ocean  west  of  Central 
America  and  Mexico,  east  of  140  degrees  West  longitude,  north  of  the  Equator  and 
generally  south  of  35  degrees  North  latitude  (OFCM  1996).  As  the  World 
Meteorological  Organization  (WMO)  Regional  Specialized  Meteorological  Center 
(RSMC)  for  Regional  Association  IV  (RA  IV),  the  National  Hurricane  Center  (NHC)  of 
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the  Tropical  Prediction  Center  (TPC)  in  Miami,  Florida,  is  responsible  for  issuing 
forecast  and  warning  advisories  for  tropical  cyclones  in  the  eastern  North  Pacific 
(Neumann  1993). 

Intensity  refers  to  the  one  minute  average  maximum  low-level  sustained  wind 
speed  of  the  tropical  cyclone  (Holland  1993).  Tropical  cyclones  in  the  eastern  North 
Pacific  are  classified  according  to  their  intensity  into  the  following  categories.  Tropical 
depressions  have  intensities  less  than  17  m  s'1,  tropical  storm  intensities  are  between  17 
m  s'1  and  32  m  s'1  and  a  hurricane’s  intensity  is  greater  than  or  equal  to  33  m  s'1  (OFCM 
1996).  With  relevant  terms  defined,  attention  can  be  turned  to  the  specific  research 
problem. 

Department  of  Defense  (DoD)  research  into  the  operational  prediction  of  tropical 
cyclone  intensity  can  be  traced  back  to  the  United  States  Navy's  Project  AROWA 
(Applied  Research;  Operational  Weather  Analyses).  Professor  Herbert  Riehl,  University 
of  Chicago,  spearheaded  research  on  this  project.  In  1951,  Project  AROWA  was  tasked 
to  develop  techniques  to  allow,  "The  prediction  of  storm  intensities  and  the  extent  of 
areas  dangerous  to  surface  ships  and  aircraft  along  the  path  of  the  tropical  cyclone"  (US 
Navy  1956).  Progress  toward  the  goal  of  predicting  tropical  cyclone  intensity  change  has 
proven  to  be  slow  and  difficult  (OFCM  1997;  AMS  1993).  Dr.  Robert  Sheets  (1990), 
former  NHC  director,  summarized  current  intensity  change  prediction  skill  as,  "sorely 
lacking."  Research  interest  in  this  topic  continues  to  the  present  and  is  the  source  of 
much  debate  (Elsberry  et  al.  1992). 
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Despite  producing  the  highest  frequency  of  tropical  cyclones  per  unit  area  on  the 
globe  (McBride  1 995),  the  eastern  North  Pacific  is  perhaps  the  least  studied  tropical 
cyclone  basin  (Tai  and  Ogura  1987).  This  unfortunate  situation  is  most  likely  due  to  a 
scarcity  of  conventional  meteorological  observations.  However,  the  lack  of  data  does  not 
diminish  a  critical  need  for  accurate  intensity  forecasts  which  impact  several  agencies. 
Only  recently  has  research  been  directed  at  understanding  the  specific  processes  affecting 
tropical  cyclone  intensity  in  the  eastern  North  Pacific  by  Whitney  and  Hobgood  (1997)  at 
The  Ohio  State  University.  Their  research  focused  on  the  relationship  between 
climatological  sea  surface  temperatures  (SST)  and  maximum  potential  intensities  (MPI) 
of  tropical  cyclones  in  this  region. 

1.2  Those  Impacted  by  the  Problem 

With  nearly  6500  km  of  coastline,  a  growing  tourism  industry  and  a  fleet  of 
fishing  vessels  threatened  by  the  devastating  impacts  of  eastern  North  Pacific  tropical 
cyclones,  the  Mexican  government  requires  accurate  intensity  forecasts  to  properly 
allocate  a  finite  amount  of  resources  to  the  associated  preparation  and  response  efforts. 
Commercial  shipping  transiting  the  region  with  destinations  west  of  the  Panama  Canal 
are  often  directly  impacted  by  tropical  cyclones. 

DoD  operations  in  the  basin  include  the  US  Navy's  Third  Fleet  based  out  of  San 
Diego,  California,  and  Pacific  Air  Forces  headquartered  in  Hawaii.  Naval  forecasters  are 
required  to  make  intensity  forecasts  out  to  72  hours  to  support  Optimum  Track  Ship 
Routing  (OTSR)  operations  (Titley  1995).  An  incorrect  intensity  forecast  can  create  a 
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potentially  dangerous  situation  for  ships.  At  a  recent  Interdepartmental  Hurricane 
Conference,  one  of  the  US  Navy’s  senior  forecasters  stressed  the  importance  of  intensity 
forecasts  to  naval  commanders  and  the  urgent  need  for  improved  accuracy  of  intensity 
forecasts  out  to  at  least  72  hours  (Barbor  1997).  Forecasters  at  the  Air  Force  Global 
Weather  Center  (AFGWC)  Americas  Region  Forecast  Branch  also  have  responsibilities 
to  support  DoD  operations  in  the  eastern  North  Pacific. 

1.3  Application  of  Satellite  Remote  Sensing  to  the  Problem 

Budgetary  constraints  prevent  the  operational  use  of  manned  aerial  reconnaissance 
to  collect  in  situ  measurements  from  tropical  cyclones  in  this  basin.  Only  on  rare 
occasions  are  research  flights  conducted  in  the  eastern  North  Pacific  (OFCM  1997). 
Therefore,  meteorological  satellites  have  become  the  main  source  of  data  for  the 
operational  specification  and  prediction  of  tropical  cyclone  intensity  (Dvorak  1990)  in 
this  otherwise  data  sparse  region.  During  a  recent  panel  discussion,  Mr.  Charles  “Chip” 
Guard,  University  of  Guam,  suggested  remotely  sensed  rainfall  measurements  might  lead 
to  advances  in  intensity  forecasting  skill  (Elsberry  et  al.  1 992). 

Three  regions  of  the  electromagnetic  spectrum  are  operationally  employed  in 
attempts  to  estimate  rainfall  from  satellites.  These  regions  are  visible  (i.e.,  X  ~  0.4  to  0.7 
pm),  infrared  (i.e.,  X  ~  10.6  to  12.6  pm)  and  microwave  (i.e.,  X  ~  0.3  to  3.0  cm).  Of 
these  spectral  regions,  only  satellite  techniques  using  microwave  measurements  can 
produce  physically  direct  estimates  of  rainfall  (Barrett  and  Beaumont  1 994).  This 
proposed  study  will  examine  the  potential  use  of  satellite  passive  microwave  rainfall 


measurements  derived  from  Special  Sensor  Microwave/Imager  (SSM/I)  radiometers 
(Hollinger  1991)  onboard  the  Defense  Meteorological  Satellite  Program  (DMSP) 
constellation  to  improve  eastern  North  Pacific  tropical  cyclone  intensity  forecasting 
techniques.  Relationships  between  parameters  obtained  from  an  operational  SSM/I-based 
rainfall  measuring  algorithm  (Ferraro  et  al.  1996)  and  12-,  24-,  36-,  48-,  60-  and  72-hour 
intensity  changes  from  postanalysis  determined  best  track  data  records  are  examined  in  an 
effort  to  identify  statistically  significant  predictors  of  intensity  change. 

Latent  heat  released  by  condensation  in  rainfall  is  recognized  as  an  important 
factor  in  tropical  cyclone  intensity  change.  Previous  studies  using  SSM/I  data  show 
correlations  between  rain  rate  or  latent  heat  release  and  intensity  trends  in  western  North 
Pacific  (Rodgers  and  Pierce  1995a;  Rao  and  MacArthur  1994)  and  Atlantic  (Rodgers  et 
al.  1994b)  tropical  cyclones.  In  this  proposed  study,  correlations  between  rainfall 
parameters  within  a  444.4  km  radius  of  the  center  and  intensity  change  are  analyzed  using 
eastern  North  Pacific  tropical  cyclone  data  from  three  years,  1992  to  1994.  Stratifications 
based  upon  tropical  cyclone  intensity,  rate  of  intensity  change,  climatology,  translation, 
landfall  and  synoptic-scale  environmental  parameters  are  studied  to  understand  factors 
that  may  affect  a  statistical  relationship  between  rainfall  parameters  and  intensity  change. 
The  predictive  skill  of  statistically  significant  rainfall  parameters  is  assessed  by  using 
independent  tropical  cyclone  data  from  another  year,  1995. 
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1.4  Goals  of  this  Research 


The  overall  goal  is  to  determine  if  SSM/I  rainfall  parameters  can  add  predictive 
skill  to  an  objective,  intensity  change  guidance  product  under  development  by  Hobgood 
(1997)  and  Petty  (1997)  at  The  Ohio  State  University.  The  resulting  guidance  product's 
skill  in  predicting  intensity  change  is  assessed  by  statistically  comparing  its  12-,  24-,  36-, 
48-,  60-  and  72-hour  intensity  change  forecast  errors  with  those  from  official  NHC, 
simple  persistence  and  no  change  forecast  techniques.  Specific  goals  related  to  the 
research  problem  are  now  defined.  Goals  listed  below  are  classified  as  primary  or 
secondary.  They  serve  as  benchmarks  for  conducting  this  research  and  are  directly  linked 
to  the  data  and  methodology  used  to  achieve  them.  Those  goals  listed  as  secondary  will 
be  fully  explored  as  time  permits. 

1.4.1  Primary  Goal  1 

Determine  if  correlations  between  rainfall  parameters  derived  from  an  operational 
algorithm  and  eastern  North  Pacific  tropical  cyclone  future  12-  to  72-hour  intensity 
changes  from  best  track  data  are  statistically  significant  at  the  0.95  level.  That  is,  a 
rainfall  parameter  will  be  considered  a  statistically  significant  predictor  of  intensity 

change  if  the  probability  that  the  correlation  coefficient  is  different  from  zero  is  greater 
than  0.95. 

1.4.2  Primary  Goal  2 

Examine  rainfall  parameter  and  intensity  change  correlations  based  upon 
stratifications  in  current  intensity,  rate  of  intensity  change,  month,  location,  SST, 
environmental  vertical  shear  of  the  horizontal  wind,  translation  speed  and  landfall. 
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1.4.3  Primary  Goal  3 

If  rainfall  parameters  are  significantly  correlated  with  intensity  change,  attempt  to 
include  them,  via  step-wise  multiple  linear  regression,  into  an  intensity  change  prediction 
model  for  eastern  North  Pacific  tropical  cyclones  already  composed  of  climatology, 
geography,  persistence  (Hobgood  1997)  and  synoptic-scale  environmental  predictors 
(Petty  1997).  Statistical  significance  level  will  be  0.95  for  entry  into  the  prediction  model 
and  0.90  to  remain. 

1.4.4  Primary  Goal  4 

If  rainfall  parameters  remain  in  the  regression-based  model,  statistically  compare 
the  12-  to  72-hour  intensity  change  prediction  errors  produced  by  the  resulting  model 
with  those  errors  produced  by  official  NHC,  Hobgood  (1997),  Petty  (1997),  simple 
persistence  and  no  change  forecast  techniques.  The  statistical  significance  of  differences 
in  errors  will  be  assessed  by  performing  paired  t  tests  at  the  0.95  level.  Success  for 
rainfall  parameters  improving  intensity  change  forecasts  is  defined  as  the  final  model 
producing  smaller  intensity  change  forecast  errors  than  the  methods  of  Hobgood  (1997) 
and  Petty  (1997).  Also,  the  final  model’s  intensity  forecast  errors  will  be  compared 
against  the  other  methods  according  to  stratifications  in  current  intensity  and  rate  of 
intensity  change. 

1.4.5  Primary  Goal  5 

Conduct  case  studies  on  individual  tropical  cyclones  to  assess  forecast  model  skill 
and  demonstrate  operational  implementation  issues.  Number  of  case  studies  conducted 
will  be  determined  by  the  time  available  for  this  research. 
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1. 4. 6  Primary  Goal  6 

Determine  if  rainfall  parameters  can  detect  the  concentric  eyewall  cycle  described 
by  Willoughby  et  al.  (1982)  and  Willoughby  (1990).  The  concentric  eyewall  cycle  is  a 
temporal  evolution  of  convective  rings  that  often  occurs  in  tropical  cyclones.  The  cycle 
has  an  important  impact  on  tropical  cyclone  intensity  change. 

1.4.7  Secondary  Goal  1 

Document  the  rainfall  characteristics  for  a  substantial  sample  of  eastern  North 
Pacific  tropical  cyclones  with  respect  to  tropical  cyclone  intensity,  rate  of  intensity 
change  and  radial  distance  from  the  center  of  circulation.  This  has  never  been  done  in  the 
eastern  North  Pacific. 

1.4.8  Secondary  Goal  2 

Analyze  the  errors  in  locating  tropical  cyclone  centers  of  circulation  with  SSM/I 
data.  Stratify  error  analysis  by  tropical  cyclone  intensity. 

1.4.9  Secondary  Goal  3 

Attempt  to  detect  a  diurnal  cycle  in  eastern  North  Pacific  tropical  cyclone  rainfall. 
Examine  the  possible  relationship  between  tropical  cyclone  intensity  and  such  a  cycle,  if 
a  cycle  is  detected. 

1.4.10  Secondary  Goal  4 

Demonstrate  the  stability  of  climatology,  geography,  persistence  and  synoptic- 
scale  environmental  variables  as  predictors  of  eastern  North  Pacific  tropical  cyclone 
intensity  change.  This  will  be  done  by  statistical  analysis  of  a  dependent  data  set  from 
different  tropical  cyclone  seasons  than  those  seasons  used  by  Hobgood  (1997)  and  Petty 


8 


(1997).  Stable  intensity  change  predictors  should  have  similar  correlation  coefficients  for 
different  tropical  cyclone  seasons. 

1.4.11  Secondary  Goal  5 

Contribute  to  the  Naval  Research  Laboratory  Marine  Meteorology  Division’s 
efforts  to  develop  a  neural  network-based  tropical  cyclone  intensity  specification 
technique  (Hawkins  et  al.  1996)  by  processing  SSM/I  data  with  adequate  coverage  of 
tropical  cyclones  in  the  eastern  North  Pacific. 

1. 4. 12  Secondary  Goal  6 

Present  the  findings  of  this  research  in  a  scholarly  journal  or  DoD  technical  report, 
at  appropriate  technical  conferences  and  in  electronic  format  on  the  Internet. 

1.5  Summary 

The  research  problem  and  associated  terms  were  defined  and  specific  goals  were 
enumerated  above.  The  main  hypothesis  underlying  this  proposed  research  is  that 
satellite  passive  microwave  measured  rainfall,  through  its  concomitant  latent  heat  release, 
will  be  highly  correlated  with  tropical  cyclone  intensity  changes  out  to  72  hours  in  the 
future.  These  strong  correlations  will  be  useful  in  improving  objective  guidance  tools  for 
the  critical  forecasting  problem  of  predicting  tropical  cyclone  intensity  changes  in  the 
eastern  North  Pacific.  Chapters  to  follow  will  provide  information  on  previous  research 
relevant  to  this  proposal,  specific  information  on  the  data  available  and  tools  to  be  used  to 
achieve  the  stated  goals.  A  final  chapter  will  discuss  accomplishments  completed  to  date 
related  to  this  proposed  research. 
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CHAPTER  2 


RELEVANT  LITERATURE  REVIEW 


2.1  Current  State  of  Tropical  Cyclone  Intensity  Forecasting 
in  the  Eastern  North  Pacific 

Improvement  of  intensity  forecasting  is  rated  as  a  “high”  priority  research 
objective  by  federal  agencies  involved  with  tropical  cyclone  research  (OFCM  1997). 
Table  1  demonstrates  the  current  intensity  forecasting  skill  for  the  eastern  North  Pacific 
during  two  time  periods  of  interest. 


Mea 

n  Absolute  Official  NHC  Intensity  Forecast  Errors  for  Eastern  North  Pacific 

Time  Period 

12  h 

ours 

24  hours 

36  hours 

48  hours 

72  hours 

m  s1 

knots 

m  s’1 

knots 

m  s‘* 

knots 

m  s"1 

knots 

m  s’1 

knots 

1988-1993* 

~ 

- 

3 

6 

- 

~ 

9 

17 

10 

20 

1990- 1 994b 

4 

7 

6 

12 

8 

16 

10 

19 

11 

22 

Table  1 :  Mean  Absolute  Official  NHC  Intensity  Forecast  Errors  (aOFCM  1 997;  bGross 
and  Lawrence  1996). 


Care  should  be  taken  when  viewing  Table  1,  because  it  demonstrates  mean  absolute 
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forecast  errors.  The  errors  for  individual  tropical  cyclones  can  be  quite  large,  especially 
for  tropical  cyclones  that  undergo  unexpected  rapid  intensity  changes.  Improvement  in 
understanding  and  predicting  intensity  change  in  the  eastern  North  Pacific  is  needed. 

2.2  Role  of  Latent  Heat  Release  in  Intensity  Change 

In  1835,  Professor  James  Pollard  Espy  proclaimed  that  the  release  of  latent  heat 
due  to  the  condensation  of  water  vapor  into  precipitation  (e.g.,  rainfall)  plays  a  major  role 
in  tropical  cyclones  (Middleton  1965).  Thus,  the  foundation  for  a  convection  or  thermal 
theory  of  tropical  cyclones  was  laid  (Khrgian  1959).  Numerical  modeling  and 
observational  studies  have  shown  that  latent  heat  release  is  the  most  important  diabatic 
process  driving  tropical  cyclones  (Anthes  1974).  Latent  heat  release  generates  available 
potential  energy  in  the  form  of  a  warm-core,  low  pressure  center  that  may  be  converted 
into  kinetic  energy.  The  conversion  to  kinetic  energy  occurs  when  friction  causes  low- 
level  winds  to  accelerate  inward  toward  the  tropical  cyclone  center.  Additionally,  this 
conversion  is  most  effective  within  about  555.5  km  of  the  center.  Therefore,  tropical 
cyclone  intensity  is  highly  sensitive  to  changes  in  latent  heat  release  (i.e.,  condensation) 
within  this  region. 

Figure  1  illustrates  the  interaction  of  the  synoptic-scale  water  and  motion  fields  in 
a  tropical  cyclone.  Moisture  (i.e.,  water  vapor)  and  liquid  water  (i.e.,  cloud  liquid  water 
and  precipitation)  and  horizontal  and  vertical  motion  are  linked  in  a  complicated, 
nonlinear  manner  via  latent  heat  release  (Atlas  and  Thiele  1981).  A  tangible  indication  of 
the  amount  of  condensation  occurring  in  a  tropical  cyclone  is  precipitation  (i.e.,  rainfall). 
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If  properly  measured,  rainfall  can  provide  useful  information  for  understanding  the  link 
between  latent  heat  release  and  tropical  cyclone  intensity  change. 


2.3  Satellite  Passive  Microwave  Measurements  of  Rainfall 

Visible  and  infrared  satellite  rainfall  estimation  methods  are  based  upon  cloud  top 
temperatures.  This  is  because  clouds,  even  very  thin  cirrus,  are  opaque  in  the  visible  and 
infrared  spectral  regions.  By  contrast,  cloud  droplets  weakly  interact  with  microwave 
electromagnetic  radiation.  However,  precipitation  strongly  affects  microwave  radiation. 
Therefore,  rainfall  can  be  readily  sensed  by  satellite-based  passive  microwave 
radiometers  in  a  physically  direct  manner.  In  crude  terms,  microwave  methods  are  able 
to  “look  through  clouds  and  see”  rainfall  occurring  below.  The  complex  interactions 
between  microwave  radiation  and  Earth’s  surface  and  atmosphere  are  summarized  in 
Figure  2.  Passive  microwave  radiometers  remotely  sense  the  microwave  radiation 
emitted,  absorbed,  reflected,  transmitted  and  scattered  by  the  Earth-atmosphere  system. 

The  physically  direct  relationship  between  microwave  radiation  and  rainfall  will 
now  be  illustrated.  The  relationship  begins  with  Planck’s  Law  for  the  radiance  or 
brightness  emitted  by  a  blackbody: 


B  (T)  = 

v 


(1) 


where  BV(T)  is  the  blackbody  emitted  radiance  or  brightness  in  units  of  W  m'2  sr'1,  h  is 
Planck’s  constant  which  is  6.6260755  *  10'34  J  s,  v  is  frequency  measured  in  units  of  Hz, 
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c  is  the  speed  of  light  in  a  vacuum  which  is  2.99792458  *  10®  m  s'1,  k  is  Boltzmann’s 
constant  which  is  1.380658  x  10'23  J  K*1  and  T is  the  absolute  temperature  with  units  of  K. 


Throughout  the  microwave  spectrum  and  troposphere,  1 ,  which  allows  the  following 


Taylor  series  approximation:  eVT'~  1  +^.  Substitution  into  Equation  1  yields: 


BIT)  - 


2  hv 


1  +-^-1 
*r 


2  ftv3  1 

c2  *1 
W 


2 v2AT  _  2kT 
c2  A2 


(2) 


The  resulting  Equation  2  is  termed  the  Rayleigh-Jeans  approximation.  The  relative  error 
of  the  approximation  is  less  than  0.002. 

Radiance  is  beneficial  for  the  remote  sensing  of  rainfall  because  it  is  conveniently 
independent  of  the  distance  to  the  rainfall,  if  the  incidence  angle  and  pathlength  are 
constant.  With  the  Rayleigh-Jeans  approximation,  radiance  is  directly  related  to  the 
absolute  or  kinetic  temperature  (i.e.,  available  from  a  thermometer)  of  a  blackbody.  Since 
space-based  radiometers  measure  radiance  instead  of  blackbody  temperature  and  a 
blackbody  is  a  theoretical  construct  not  experienced  in  nature,  it  is  useful  to  compare  the 
radiance  of  a  real  emitter  or  graybody  to  that  of  a  blackbody: 

2*r/ 

l a  =  _  £*  =  e 

Bv(T )  2kT  T  T  ~  Ey  ^ 

A2 

The  ratio  in  Equation  3  quantifies  how  efficiently  a  graybody  emits  microwave  radiation 
as  compared  to  a  blackbody  with  the  same  absolute  temperature.  B  \(T)  is  the  radiance  of 
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a  graybody  which  is  remotely  sensed.  iT'is  the  absolute  temperature  a  blackbody  would 
have  for  a  radiance  of  B  \(T).  T  is  again  the  absolute  temperature  of  the  graybody.  TB  is 
defined  as  the  brightness  temperature  and  is  related  to  the  absolute  temperature,  T,  by  a 
physical  property  of  the  graybody  called  emissivity  or  emittance,  ev. 


transmitted  radiance  at  v 
incident  radiance  at  v 


(4) 


Pv 


reflected  radiance  at  v 
incident  radiance  at  v 


(5) 


absorbed  radiance  at  v 
incident  radiance  at  v 


(6) 


Other  useful  properties  in  the  discussion  of  microwave  radiance  are  defined  in 
Equations  4  to  6  as  transmittance  (i.e.,  tv),  reflectance  (i.e.,  pv)  and  absorptance  (i.e.,  av). 
These  are  related  by  the  following  identity  due  to  conservation  of  energy  considerations: 

"*■  Pv +  =  1  •  Kirchhoff  s  Law  states  that  in  the  case  of  local  thermodynamic 

equilibrium,  a  valid  assumption  for  the  Earth’s  atmosphere  below  100  km,  av  =  ev. 

An  equation  describing  the  transfers  of  microwave  radiation  illustrated  in  Figure  2 
can  be  written  in  terms  of  TB  and  T.  Using  the  plane-parallel  geometry  shown  in  Figure  3 
and  neglecting  scattering,  an  equation  for  microwave  radiative  transfer  involving  rainfall 
is  derived  (Grody  1997;  Kidder  and  Vonder  Haar  1995): 

dT 

COS0T£  =  a°(r-r*>  (7) 
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The  volume  absorption  coefficient,  oa,  has  units  of  m'1.  Above  the  rain  layer,  o2  =  0  m'1. 
Therefore,  integration  of  Equation  7  is  simplified  to  the  depth  of  the  rain  layer,  D.  Also, 
the  integration  can  be  further  simplified  by  assuming  that  the  temperature  of  the  rain 
layer,  Tn  is  constant.  The  incidence  angle,  0,  is  measured  from  local  vertical  to  the 
radiometer.  Equation  8  below  shows  that  TB  measured  at  the  radiometer  can  be 
represented  as  the  sum  of  the  four  components  shown  graphically  in  Figure  3: 

tb  a  Ti  +  T2  +  T3  +  (8) 

Each  term  on  the  right  side  of  Equation  8  is  defined  below: 

Tx  =  r  etsec6  (9) 


T2  =  (1  -t)Tc  T2sec6 


(10) 


T '  =  (1 


•e)7r(l  -Tsec0)Tsec6 


74  =  Tr(  1  -t“c0) 


(11) 

(12) 


T,  represents  the  microwave  emission  by  the  surface  with  temperature  Ts  and 
transmission  through  the  atmosphere.  Cosmic  microwave  transmission  through  the 
atmosphere  and  reflection  at  the  surface  is  described  by  T2.  Tc  is  the  cosmic  radiometric 
temperature  which  is  quite  small  (i.e.,  Tc  =  2.7  K).  The  rain  layer’s  downwelling 
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microwave  emission,  transmission  through  the  atmosphere  and  reflection  by  the  surface  is 
given  by  T3.  T4  describes  the  upwelling  microwave  emission  by  the  rain  layer  and 
transmission  through  the  atmosphere.  Transmittance  of  the  rain  layer  is  approximated  by: 
t  =  e (  °  \  Substitution  of  Equations  9  to  12  into  Equation  8  and  rearrangement  of  terms 
yields: 

Tb  *  r  exsec6+(l  -8)rcT2sec6  +  rr[l  -exsec0-(i  -£)x2sec6]  (13) 

If  the  small  cosmic  term  is  neglected  in  Equation  13,  the  impact  of  rainfall  on  TB 
can  be  simply  illustrated.  For  a  nonraining  atmosphere,  x  ~  1,  which  means  TB  ~  Ts  e.  As 
rain  rate  increases,  oa  increases,  x  approaches  0  and  TB  converges  to  Tr.  Therefore,  Tg 
measurements  can  provide  extremely  useful  information  by  their  response  to  rainfall. 

Some  problems  confound  the  otherwise  simple  relationship  between  TB  and 
rainfall.  Over  water,  TB  increases  rapidly  with  rain  rate  because  s  values  for  water  are 
relatively  low  (i.e.,  e  «  0.5).  Thus,  water  provides  a  radiometrically  cool  background 
compared  to  a  warm  rain  layer  target.  However,  £  values  over  land  are  relatively  high 
(i.e.,  e  *  0.9).  As  a  result,  land  is  a  warm  background  for  a  warm  rain  layer.  Therefore,  a 
change  in  TB  due  to  rainfall  emissions  over  land  is  usually  too  small  to  detect.  Another 
problem  is  that  TB  depends  on  the  product  of  oa£>.  However,  rain  rate  is  highly  related  to 
oa  and  independent  of  D.  Therefore,  the  depth  of  the  rain  layer,  D,  needs  to  be  specified 
in  order  for  TB  measurements  to  be  used  for  explicitly  determining  rain  rate.  Also,  oa 
depends  not  only  on  rain  rate,  but  cloud  liquid  water  and  water  vapor,  as  well. 

Consequently,  cloud  liquid  water  and  water  vapor  information  are  necessary  to  perform 
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calculations  of  rain  rate  from  measurements  of  TB. 

By  neglecting  scattering  and  cosmic  microwave  radiation,  Equation  13  implies 
that  as  rain  rate  increases,  TB  increases  to  a  limit  of  Tr.  However,  in  reality  there  is  a  point 
at  which  increasing  rain  rate  actually  results  in  decreasing  values  of  TB.  This  is  because 
as  rain  rate  increases,  the  size  of  hydrometeors  in  the  rain  layer  (e.g.,  very  large  rain 
drops,  ice,  graupel  and  hail)  increases  and  approaches  the  wavelength  of  microwave 
radiation.  Scattering  then  occurs  which  obscures  emission  by  the  rain  layer  and  surface 
and  allows  reflection  of  extremely  cold  cosmic  microwave  radiation  off  the  top  of  the  rain 
layer  into  the  radiometer’s  field  of  view.  The  resulting  reduction  in  TB  occurs  over  water 
and  land  surfaces. 

The  previous  discussion  in  this  section  describes  two  regimes  for  satellite  passive 
microwave  measurements  of  rainfall:  emission  and  scattering.  Up  to  a  threshold  rain  rate 
governed  by  wavelength,  emission  by  rainfall  over  water  allows  TB  measurements  to  be 
used  to  calculate  rain  rate.  Beyond  that  threshold  rain  rate,  scattering  by  large 
hydrometeors  permits  TB  measurements  to  be  used  to  infer  rain  rate  over  water  and  land. 
In  general,  microwave  frequencies  below  20  GHz  are  used  in  emission-based  algorithms. 
At  frequencies  above  60  GHz,  scattering  effects  are  utilized  by  passive  microwave 
algorithms.  Between  these  frequencies,  emission  and  scattering  effects  are  found.  An 
excellent  summary  of  the  current  state-of-the-art  algorithms  developed  by  researchers  to 
calculate  rain  rate  from  passive  microwave  brightness  temperatures  was  written  by 
Wilheit  et  al.  (1994).  The  algorithm  to  be  used  in  this  proposed  research  utilizes  both 
emission  and  scattering  regimes  to  measure  rainfall  in  eastern  North  Pacific  tropical 
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cyclones.  Details  of  the  relevant  algorithm  and  reasons  it  was  selected  for  this  research 
are  discussed  below. 

2.4  Rainfall  Algorithm  to  be  Used  in  this  Research 

An  important  distinction  in  rainfall  types  found  in  tropical  cyclones  is  convective 
versus  stratiform  (Houze  1993).  Convective  precipitation  occurs  in  cumulus  and 
cumulonimbus  clouds  of  the  eyewall  and  convective  rain  bands.  Nimbostratus  clouds  in 
tropical  cyclones  produce  stratiform  rainfall.  Tropical  rain  rates  less  than  or  equal  to  1 
mm  hr'1  are  mostly  stratiform  and  those  greater  than  10  mm  hr'1  are  generally  convective 
(Tokay  and  Short  1996).  Rain  rates  between  these  two  values  can  be  caused  by  both 
types  of  rainfall.  However,  a  threshold  between  the  two  types  at  6  mm  hr'1  was  proposed 
by  Johnson  and  Hamilton  (1988). 

Processes  governing  convective  precipitation  take  place  above  the  freezing  level 
and  are  termed  “cold  cloud.”  Passive  microwave  algorithms  best  suited  for  this  type  of 
precipitation  measurement  are  those  that  take  advantage  of  the  pronounced  scattering 
effects  of  the  resulting  precipitation-sized  ice  particles.  The  “warm  cloud”  processes  of 
stratiform  rainfall  take  place  at  or  below  the  freezing  level.  Scattering-based  algorithms 
are  not  useful  in  this  warm  cloud  regime.  However,  emission-based  algorithms  can 
detect  stratiform  rainfall  over  water. 

Greater  than  90  percent  of  the  rainfall  associated  with  a  tropical  cyclone  occurs 
within  222.2  km  of  the  center  of  circulation  (Hughes  1952;  Riehl  1954).  Of  this  rainfall, 
approximately  90  percent  of  the  areal  coverage  within  this  region  is  stratiform  (Jorgensen 
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1984).  Due  to  higher  rain  rates,  convective  rainfall  contributes  about  60  percent  of  the 
total  rainfall  (Marks  1990;  Jorgensen  1984).  Therefore,  the  measurement  of  both 
convective  and  stratiform  rainfall  is  important  to  this  study. 

The  rainfall  algorithm  chosen  for  this  proposed  research  was  developed  by  the 
National  Environmental  Satellite,  Data  and  Information  System  (NESDIS)  Office  of 
Research  Applications  (ORA).  It  was  developed  to  determine  instantaneous  rain  rates 
using  brightness  temperature,  TB,  data  from  the  DMSP  SSM/I  (Ferraro  et  a!.  1996). 
Specific  details  concerning  DMSP  SSM/I  data  will  be  provided  in  Chapter  3.  There  are 
two  main  components  of  the  NESDIS/ORA  over  water  algorithm:  scattering  and 
emission.  The  scattering  component  relates  a  scattering  index  (Grody  1991;  Ferraro  et  al. 
1 994)  to  rain  rates.  This  scattering  index  is  the  difference  between  the  actual  vertically 
polarized  brightness  temperature  at  85.5  GHz  and  one  predicted  for  a  scattering-free 
atmosphere  by  using  the  vertically  polarized  brightness  temperatures  at  19.4  and  22.2 
GHz.  A  brightness  temperature  screening  process  is  used  to  prevent  false  rain  retrievals 
(e.g.,  sea  ice).  The  emission  component  is  calculated  by  using  cloud  liquid  water  values 
greater  than  0.20  mm  estimated  using  the  vertically  polarized  brightness  temperatures  at 
19.4,  22.2  and  37.0  GHz  (Weng  and  Grody  1994).  Scattering  index  and  cloud  liquid 
water  relationships  with  rain  rate  were  empirically  developed  using  coincidental  ground- 
based  radar  rainfall  measurements  from  around  the  world  (Ferraro  and  Marks  1995).  The 
range  of  rain  rates  detected  by  the  algorithm  is  between  0.2  mm  hr'1  and  35.0  mm  hr'. 
Validation  of  the  algorithm  found  a  relative  error  between  SSM/I  and  radar  rain  rates  of 
0.10  to  0.20.  Finally,  a  practical  reason  for  using  the  NESDIS/ORA  algorithm  is  that  it 
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was  selected  by  the  Shared  Processing  Program  Algorithm  Research  Panel  as  the  SSM/I 
rainfall  algorithm  to  be  operationally  implemented  by  DoD  and  other  federal  agencies 
(Colton  and  Poe  1994;  Ferraro  et  al.  1996).  Therefore,  this  rainfall  product  is  readily 
available  to  NHC  forecasters  via  the  Shared  Processing  Program  and  the  Internet. 

2.5  Satellite  Passive  Microwave  Measurements  Related  to 
Tropical  Cyclone  Intensity 

Links  between  space-based  passive  microwave  measurements  related  to 
precipitation  and  tropical  cyclone  intensity  have  been  explored  for  over  20  years.  The 
first  research  on  this  specific  topic  was  conducted  by  Adler  and  Rodgers  (1977).  By 
examining  six  observations  of  western  North  Pacific  tropical  cyclone  Nora  (1972),  they 
found  latent  heat  release  (LHR)  increased  as  intensity  increased.  LHR  (in  units  of  W) 
within  a  444.4  km  radius  of  the  tropical  cyclone  center  was  calculated  with  a  rainfall 
algorithm  using  the  19.4  GHz  horizontally  polarized  brightness  temperatures  measured 
by  the  Nimbus-5  Electrically  Scanning  Microwave  Radiometer  (ESMR-5). 

Rodgers  and  Adler  (1981)  later  used  71  observations  by  the  ESMR-5  of  21 
tropical  cyclones  during  the  1973  to  1975  seasons.  49  of  the  observations  were  for  18 
western  North  Pacific  tropical  cyclones  and  22  of  the  observations  were  for  3  eastern 
North  Pacific  tropical  cyclones.  A  correlation  of  0.71,  significant  at  the  0.99  level,  was 
found  between  LHR  within  a  444.4  km  radius  and  current  intensity.  One  case  study  of 
eastern  North  Pacific  tropical  cyclone  Doreen  (1973)  demonstrated  their  finding  that 
maximum  LHR  generally  occurred  one  to  two  days  prior  to  maximum  intensity.  This 
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finding  confirmed  earlier  numerical  tropical  cyclone  prediction  model  results  that  showed 
maximum  intensity  occurred  one  to  three  days  after  maximum  LHR  (Kurihara  and  Tuleya 
1974;  Rosenthal  1978).  Of  further  interest,  Rodgers  and  Adler  (1981)  found  that  the 
three  eastern  North  Pacific  tropical  cyclones  were  more  compact  and  had  less  rainfall 
than  those  of  the  western  North  Pacific. 

Encouraged  by  the  results  of  Adler  and  Rodgers  (1977)  and  Rodgers  and  Adler 
(1981),  Hunter  et  al.  (1981)  made  the  first  attempt  at  an  objective  tropical  cyclone 
intensity  forecast  guidance  tool  including  microwave  brightness  temperature  data.  They 
employed  an  eigenvector  representation  of  the  ESMR-5  19.4  GHz  horizontally  polarized 
brightness  temperature  images  centered  on  120  observations  of  29  tropical  cyclones  from 
the  1973  to  1974  seasons.  25  were  tropical  cyclones  in  the  western  North  Pacific  and  4 
occurred  in  the  eastern  North  Pacific.  By  adding  ESMR-5  information  to  climatology, 
geography  and  persistence  information  from  the  best  track  data,  they  were  able  to 
increase  the  explained  variance  of  intensity  by  0.07  at  the  24-hour  forecast  period  and  by 
0. 16  at  the  72-hour  forecast  period,  when  compared  to  using  only  best  track  data.  Then- 
guidance  product  based  upon  best  track  and  ESMR-5  predictors  produced  mean  absolute 
forecast  errors  for  the  24-  and  72-hour  periods  that  were,  respectively,  2  and  4  m  s'1  less 
than  the  mean  absolute  official  forecast  errors  by  the  Joint  Typhoon  Warning  Center 
(JTWC)  for  the  same  forecast  periods  and  tropical  cyclone  seasons. 

With  the  launch  of  the  first  SSM/I  in  1987,  routine  operational  microwave 
measurements  of  tropical  cyclones  from  satellites  became  possible.  The  first  researcher 
to  relate  microwave  measurements  from  the  SSM/I  to  tropical  cyclone  intensity  was 
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Rhudy  (1989).  He  used  23  observations  of  12  tropical  cyclones  in  the  western  North 
Pacific  and  Atlantic  during  1987  to  1988.  LHR  within  a  444.4  km  radius  and  current 
intensity  had  a  correlation  of  0.62,  significant  at  the  0.99  level.  85.5  GHz  vertically 
polarized  brightness  temperature  differences  within  a  333.3  km  radius  yielded 
correlations  of  0.64  and  0.74,  significant  at  the  0.99  level,  with  current  and  future  24- 
hour  intensities,  respectively. 

Relationships  between  SSM/I-measured  rainfall  and  future  24-hour  intensity 
change  were  examined  for  27  observations  of  12  tropical  cyclones  in  the  western  North 
Pacific  and  Atlantic  during  1987  to  1988  (MacArthur  1991;  Rao  etal.  1991;  Rao  and 
MacArthur  1994).  They  found  a  0.68  correlation,  significant  at  the  0.95  level,  between 
average  volumetric  rain  rate  (in  units  of  mm3  hr'1)  within  a  222.2  km  radius  and  future  24- 
hour  intensity  change.  Case  studies  of  western  North  Pacific  tropical  cyclones  Lynn 
(1987)  and  Dinah  (1987)  yielded  extremely  high  correlations  between  average  volumetric 
rain  rate  within  a  222.2  km  radius  and  future  24-hour  intensity  change.  The  correlation 
for  Lynn’s  seven  observations  was  0.97  and  Dinah’s  four  observations  was  0.99.  Overall, 
they  found  a  differential  effect  of  rainfall  on  future  intensity  change  due  to  location  of 
rainfall  with  respect  to  the  center  of  circulation.  Within  a  222.2  km  radius  from  the 
center,  rainfall  had  a  highly  positive  correlation  with  future  24-hour  intensity  change. 
However,  within  the  annular  region  between  222.2  to  444.4  km  from  the  center,  rainfall 
had  a  weakly  negative  correlation  with  future  24-hour  intensity  change.  It  was  noted  that 
if  a  rain  rate  maximum  occurred  outside  of  the  core  region  (i.e.,  beyond  a  222.2  km 
radius),  short-term  weakening  of  intensity  could  be  expected. 
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McCoy  (1991)  used  85.5  GHz  vertically  polarized  brightness  temperature 
differences  between  regions  (i.e.,  outer-versus-inner  and  lefit-versus-right)  within  a  444.4 
km  radius  to  find  a  0.61  correlation  with  future  24-hour  intensity  change.  He  also  noted 
that  the  appearance  of  convective  rain  bands  beyond  222.2  km  from  the  center  usually 
preceded  a  weakening  of  future  24-hour  intensity.  His  study  was  based  on  25 
observations  of  western  North  Pacific  and  Atlantic  tropical  cyclones  during  1987  to  1988. 
In  this  study,  McCoy  (1991)  was  the  first  researcher  to  suggest  that  SSM/I  measurements 
might  be  able  to  monitor  the  convective  ring  or  concentric  eyewall  cycle  discussed  by 
Willoughby  et  al.  (1982  and  1984)  and  Willoughby  (1988  and  1990).  This  convective 
ring  cycle  is  a  spatial  and  temporal  evolution  of  axisymmetric  convective  rainfall  rings 
that  commonly  occurs  in  symmetric  tropical  cyclones  with  intensities  greater  than  or 
equal  to  45  m  s'1.  Often  a  single  vigorous  circular  ring  of  rainfall  or  an  eyewall  within  a 
1 1 1 . 1  km  radius  will  contract  inward  with  a  concomitant  intensity  increase.  Sometimes 
concentric  outer  convective  rings  will  form  between  a  1 1 1.1  to  222.2  km  radius  and 
contract  around  the  existing  eyewall.  This  outer  convective  ring  causes  the  inner  eyewall 
to  dissipate  and  intensity  decreases.  As  the  outer  convective  ring  continues  to  contract  it 
often  supplants  the  original  eyewall  and  intensity  increases  when  this  new  eyewall 
contracts.  The  entire  eyewall  cycle  normally  takes  between  12  to  36  hours  to  complete. 
Multiple  concentric  eyewall  cycles  may  occur  provided  the  tropical  cyclone  does  not 

make  landfall,  experience  excessive  environmental  wind  shear  or  travel  over  relatively 
cold  ocean. 
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Felde  and  Glass  (1991)  used  the  fractional  coverage  of  85.5  GHz  horizontally 
polarized  brightness  temperatures  less  than  220  K  within  a  tropical  cyclone-centered 
circular  area  with  a  55.5  km  radius.  Their  study  used  17  observations  of  1 1  western 
North  Pacific  tropical  cyclones  during  the  1987  season.  The  correlation  between 
fractional  coverage  of  intense  convection,  indicated  by  depressed  85.5  GHz  brightness 
temperatures,  and  current  intensity  was  0.73.  Also,  the  difference  between  the  19.4  GHz 
horizontally  and  vertically  polarized  brightness  temperatures  (i.e.,  an  indication  of 
microwave  emission  by  rainfall)  within  the  same  55.5  km  radius  yielded  a  0.77 
correlation  with  current  intensity.  For  nine  observations  of  intensifying  tropical  cyclones 
the  correlation  increased  to  0.85  and  for  eight  observations  of  weakening  tropical 
cyclones  the  correlation  was  0.70.  A  later  study  by  Glass  and  Felde  (1992)  used  25 
observations  of  19  western  North  Pacific  tropical  cyclones  during  1987  to  1988.  The 
same  85.5  GHz  parameter  described  above  produced  a  0.93  correlation  with  current 
intensity  for  15  intensifying  tropical  cyclones  and  a  0.70  correlation  for  10  weakening 
tropical  cyclones. 

Alliss  et  al.  (1992)  conducted  a  case  study  of  Atlantic  tropical  cyclone  Hugo 
(1989)  using  nine  SSM/I  observations.  They  found  SSM/I-derived  rainfall  parameters 
calculated  for  areas  within  a  444.4  km  radius  were  highly  correlated  with  current 
intensity.  The  highest  correlation  coefficient  was  0.93,  significant  at  the  0.99  level.  It 
was  suggested  that  Hugo’s  intensity  changes  might  be  related  to  the  convective  ring 
cycle.  In  a  similar  case  study  of  Atlantic  tropical  cyclone  Florence  (1988)  using  five 
observations,  Alliss  et  al.  (1993)  found  the  average  rain  rate  (in  units  of  mm  hr*1)  within  a 
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1 1 1 . 1  km  radius  had  a  correlation  of  0.93  with  current  intensity,  significant  at  the  0.98 
level.  Both  of  these  studies  found  that  as  the  azimuthally  averaged  rainfall  maximum 
moved  inward,  intensity  increased. 

In  a  previous  study,  this  author  used  17  SSM/I  observations  of  western  North 
Pacific  tropical  cyclone  Freda  (1987),  climatological  sea  surface  temperature  and  global 
numerical  weather  prediction  (NWP)  model  data  to  understand  the  relationships  between 
rainfall,  environmental  forcing  and  intensity  (West  1993a  and  b).  Changes  in  LHR 
preceded  changes  of  intensity  by  about  12  hours.  Moisture  flux  convergence  in  the 
middle  troposphere  initiated  convective  rings  in  the  outer  core  (i.e.,  within  1 1 1.1  to  222.2 
km).  Interactions  with  upper-level  troughs  aided  inner  core  (i.e.,  within  111.1  km) 
rainfall,  if  vertical  wind  shear  was  weak.  Changes  in  intensity  were  related  to  the 
convective  ring  cycle. 

A  case  study  of  western  North  Pacific  tropical  cyclone  Flo  (1990)  using  seven 
SSM/I  observations  found  a  correlation  of  0.70  between  rainfall  within  a  111.1  km  radius 
and  current  intensity  (Melton  1994).  Zhao  (1994)  used  three  cases  of  western  North 
Pacific  tropical  cyclones  from  1992  to  relate  rainfall  parameters  to  future  intensity.  She 
studied  a  total  of  37  SSM/I  observations  from  21,  7  and  9  observations  of  tropical 
cyclones  Gay,  Elsie  and  Hunt,  respectively.  For  intensifying  tropical  cyclones,  the 
correlation  between  a  measure  of  total  (i.e.,  convective  and  stratiform)  rainfall  within  a 
111.1  km  radius  and  future  12-hour  intensity  change  was  0.71,  highly  significant  at  the 
0.999  level.  A  convective  rainfall  parameter  within  a  222.2  km  radius  yielded  a  0.63 
correlation,  significant  at  the  0.91  level,  with  future  12-hour  intensity  change  in 
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weakening  tropical  cyclones.  Zhao  (1994)  also  found  the  convective  ring  cycle  occurred 
in  these  cases  and  maximum  rainfall  occurred  prior  to  maximum  intensity. 

Rodgers  et  al.  (1994b)  studied  18  North  Atlantic  tropical  cyclones  from  1987  to 
1989  with  103  SSM/I  observations.  They  found  that  more  intense  tropical  cyclones  had 
more  rainfall.  Additionally,  correlations  between  past  12-hour  trends  in  both  LHR  within 
a  1 1 1 . 1  km  radius  and  intensity  varied  according  to  the  current  intensity.  That  is,  for 
tropical  depressions  the  correlation  was  0.06,  for  tropical  storms  it  was  -0.27  and  for 
hurricanes  it  was  0.78,  all  significant  at  the  0.99  level.  Rodgers  and  Pierce  (1995a)  later 
examined  western  North  Pacific  tropical  cyclones  from  1987  to  1992  with  257  SSM/I 
observations.  For  123  tropical  depression  observations  the  correlation  between  past  12- 
hour  trends  in  both  average  rain  rate  within  a  1 1 1.1  km  radius  and  intensity  was  0.25,  for 
61  tropical  storm  observations  it  was  0.51  and  for  73  typhoon  observations  it  was  0.04, 
all  significant  at  the  0.99  level.  Additionally,  a  diurnal  variation  of  rainfall  within  a  444.4 
km  radius  was  found  for  weaker  tropical  cyclones  (i.e.,  less  than  tropical  storm  strength). 
An  early  morning  maximum  and  an  evening  minimum  were  found  for  these  weaker 
tropical  cyclones.  For  stronger  tropical  cyclones,  especially  within  the  inner  core,  no 
significant  diurnal  variation  was  found.  These  findings  support  the  results  of  a  study  by 
Hobgood  (1986)  that  included  the  diurnal  cycle  of  net  radiation  at  the  cloud  tops  in  a 
numerical  tropical  cyclone  prediction  model.  Rodgers  et  al.  (1994b)  and  Rodgers  and 
Pierce  (1995a)  both  found  that  higher  rain  rates  occurred  near  the  center  of  tropical 
cyclones  and  that  the  spatial  and  temporal  change  of  azimuthally  averaged  rainfall 
demonstrated  the  convective  ring  cycle. 
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Using  SSM/I,  climatological  SST  and  global  NWP  model  data,  Rodgers  et  al. 
(1994a)  and  Rodgers  and  Pierce  (1995b)  examined  case  studies  of  tropical  cyclones  to 
understand  environmental  forcing,  rainfall  and  intensity  relationships.  Rodgers  et  al. 
(1994a)  used  30  SSM/I  observations  of  3  Atlantic  tropical  cyclones  during  1989.  Tropical 
cyclone  Dean  (1989)  had  7  observations,  Gabrielle  had  1 1  and  Hugo  had  12.  In  Rodgers 
and  Pierce  (1995),  18  observations  of  western  North  Pacific  tropical  cyclone  Bobbie 
(1992)  were  used.  In  both  studies  and  all  cases,  the  convective  ring  cycle  played  an 
important  role  in  intensity  change.  A  convective  ring  often  formed  in  the  outer  core  (i.e., 
between  111.1  and  222.2  km  radii)  and  propagated  inward.  Maximum  intensity  occurred 
when  the  ring  reached  the  inner  core  (i.e.,  within  1 1 1 . 1  km  of  the  center).  The  formation 
of  outer  convective  rings  usually  led  to  short-term  weakening.  Environmental  forcing 
affects  were  similar  to  those  found  by  West  (1993a  and  b).  In  a  recent  case  study  by 
Rodgers  et  al.  (1997),  13  SSM/I  observations  of  LHR  within  a  55.5  km  radius  of  Atlantic 
tropical  cyclone  Opal  (1995)  were  studied  to  understand  the  rapid  intensity  change  that 
occurred.  Three  convective  bursts  took  place  prior  to  a  rapid  intensity  increase.  They 
proposed  that  the  rapid  intensity  change  was  due  to  the  massive  amounts  of  latent  heat 
released  by  the  convective  bursts. 

Recently,  Cecil  (1997)  related  85.5  GHz  brightness  temperature-based  ice 
scattering  signatures  within  a  1 1 1.1  km  radius  to  future  24-hour  intensity  for  90  SSM/I 
observations  of  western  and  eastern  North  Pacific  and  Atlantic  tropical  cyclones.  The 
correlation  was  0.68  for  all  basins.  For  the  17  observations  of  6  eastern  North  Pacific 
tropical  cyclones  the  correlation  improved  to  0.89. 
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Previous  studies  of  satellite  microwave  measurements  of  tropical  cyclones 
indicate  that  parameters  related  to  rainfall  are  highly  correlated  with  tropical  cyclone 
intensity  change.  Used  in  concert  with  other  factors  related  to  tropical  cyclone  intensity 
change,  SSM/I-derived  rainfall  parameters  ought  to  aid  in  the  prediction  of  eastern  North 
Pacific  tropical  cyclone  intensity  change. 

2.6  Factors  Related  to  Tropical  Cyclone  Intensity  Change 

Strong  correlations  between  rainfall  parameters  and  tropical  cyclone  intensity 
were  shown  above.  However,  there  are  other  factors  affecting  tropical  cyclone  intensity 
change  that  might  confound  attempts  to  develop  a  predictive  tool  for  intensity  change 
based  solely  on  rainfall  parameters.  A  brief  summary  of  research  on  these  factors  is 

reviewed  below.  These  factors  will  be  used  in  the  stratified  statistical  analysis  discussed 
in  Section  3.10. 

The  current  intensity  of  a  tropical  cyclone  can  affect  the  relationship  between 
LHR  and  future  intensity  change.  Studies  by  Shapiro  and  Willoughby  (1982),  Rodgers  et 
al.  (1994b)  and  Rodgers  and  Pierce  (1995a)  showed  that  as  tropical  cyclone  intensity 
increased,  the  correlation  between  LHR  within  the  core  region  and  future  intensity  change 
became  stronger.  This  increasingly  efficient  relationship  between  LHR  and  intensity 

change  for  more  intense  tropical  cyclones  was  attributed  to  increasing  lower-tropospheric 
inertial  stability. 

Rate  of  intensity  change  is  an  important  factor  in  understanding  LHR  and  future 
intensity  change.  Rapid  intensity  changes  (i.e.,  greater  than  11  ms'1  over  a  12-hour 
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period)  are  extremely  hard  to  forecast  and  cause  large  official  intensity  change  forecast 
errors  (Elsberry  et  al.  1992;  Mundell  1990  and  1991;  Sampson  et  al.  1995).  Mundell 
(1990  and  1991)  found  that  extremely  intense  convection,  inferred  by  infrared  satellite 
measurements,  within  a  666.6  km  radius  of  the  tropical  cyclone  center  usually  preceded 
rapid  intensification  events  by  12  hours. 

The  climatological  (i.e.,  geographical  and  monthly)  variability  of  tropical  cyclone 
intensity  characteristics  was  first  studied  by  Frank  and  Jordan  (1960).  Whitney  (1995) 
and  Whitney  and  Hobgood  (1997)  conducted  such  a  study  using  3 1  years  of  data  for  the 
eastern  North  Pacific.  In  this  study,  geographical  stratifications  were  north  or  south  of  1 8 
degrees  North  latitude  and  east  or  west  of  1 1 0  degrees  West  longitude.  Significant 
geographical  and  monthly  variations  of  intensity  were  found. 

A  necessary ,  but  not  sufficient,  condition  for  tropical  cyclone  formation  is  a  sea 
surface  temperature  of  at  least  26  °C  (Palmen  1948).  Also,  in  a  study  of  eastern  North 
Pacific  tropical  cyclones  by  Whitney  (1995)  and  Whitney  and  Hobgood  (1997),  the 
maximum  potential  intensity  was  directly  related  to  climatological  SST. 

Weak  tropospheric  vertical  shear  of  the  horizontal  wind  is  another  necessary,  but 
not  sufficient,  condition  for  tropical  cyclone  development  (Gray  1968).  The  vector 
difference  of  the  horizontal  wind  at  200  and  850  mb  is  often  used  to  quantify  this  vertical 
shear.  In  a  study  of  western  North  Pacific  tropical  cyclones,  Zehr  (1992)  found  values  of 
200  to  850  mb  vertical  shear  greater  than  15  ms'1  to  be  excessive  and  not  conducive  to 
tropical  cyclone  development. 
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Translation  speeds  greater  than  10  m  s'1  are  not  conducive  to  intensification 
(Sampson  et  al.  1995).  Tropical  cyclone  landfall  has  complex  interactions  with  intensity 
change  (Merrill  1987).  Therefore,  this  proposed  research  will  exclude  data  at  landfall. 

2.7  Statistical  Prediction  of  Eastern  North  Pacific  Tropical  Cyclone 
Intensity  Change 

Recently,  research  at  The  Ohio  State  University  has  been  directed  at  developing 
statistically  based  guidance  for  forecasting  intensity  change  in  eastern  North  Pacific 
tropical  cyclones  (Hobgood  1997;  Petty  1997).  The  research  by  Hobgood  (1997) 
parallels  and  expands  upon  similar  research  in  the  North  Atlantic  by  Jarvinen  and 
Neumann  (1979).  Using  ten  years  (i.e.,  1982  to  1993)  of  best  track  data  (see  Section  3.1 
for  details  on  this  type  of  data)  and  climatological  SST,  ten  predictors  of  future  intensity 
change  out  to  72  hours  were  identified.  These  forecast  variables  are  summarised  in  Table 
2.  Also  shown  in  Table  2  is  the  variance  of  intensity  change  explained  by  these  variables 
at  each  forecast  period  and  the  mean  absolute  forecast  errors  for  a  regression  model  tested 
against  independent  best  track  data  from  1 994.  These  errors  compare  favorably  with 
mean  absolute  official  NHC  intensity  forecast  errors. 

Petty  s  (1997)  research  adds  synoptic-scale  environmental  forcing  information 
derived  from  global  NWP  model  data  (see  Section  3.1  for  more  details)  to  the  statistical 
guidance  product  developed  by  Hobgood  (1997).  This  research  is  analogous  to  and  goes 
beyond  that  of  DeMaria  and  Kaplan’s  (1994)  study  in  the  North  Atlantic.  Additionally, 
Petty  (1997)  uses  real-time,  instead  of  climatological,  SST.  By  combining  the  predictors 
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derived  from  the  NWP  model  and  best  track  data  for  three  years  (i.e.,  1989  to  1991),  he 
identified  a  total  of  13  predictors  of  intensity  change.  His  wind  shear,  thermal,  moisture 
and  eddy  flux  forecast  variables  are  summarized  in  Table  2.  Also  shown  in  Table  2  is  the 
variance  of  intensity  change  explained  by  these  variables  at  each  forecast  period. 


Summary  of  Eastern  North  Pacific  Tropical  Cyclone  Intensity  Change 
Forecasting  Techniques  Under  Development  at  The  Ohio  State  University 


Climatology,  geography  and  persistence  forecast  variables  identified  by  Hobgood  (1997)  in 
decreasing  order  of  importance: 

1.  Latitude*;  2.  Current  intensity*;  3.  Previous  12  hour  intensity  trend*;  4.  Longitude;  5.  Distance  to  land 
calculated  with  Kaplan  (1992);  6.  Absolute  value  of  the  Julian  date  minus  237*;  7.  Potential 
intensification  calculated  by  subtracting  the  current  intensity  from  the  maximum  potential  intensity 
defined  by  Whitney  and  Hobgood  (1997);  8.  Zonal  component  of  translation;  9.  Meridional  component 
of  translation;  10.  Translation  speed*. 


Forecast  Period 


12  hour 


24  hour 


36  hour 


48  hour 


60  hour 


72  hour 


Explained  Variance  of  Intensity  Change  (r2)  for  Hobgood  (1997) 


1 994  Mean  Absolute  Forecast  Errors  (m  s  ')  for  Hobgood  (1997) 


Synoptic-scale  wind  shear,  thermal,  moisture  and  eddy  flux  forecast  variables  identified  by  Petty 
(1997)  in  decreasing  order  of  importance: 

1.  Absolute  vertical  wind  shear  between  500  and  850  mb;  2.  Previous  24  hour  trend  of  absolute  vertical 
wind  shear  between  200  and  850  mb;  3.  Zonal  component  of  vertical  wind  shear  between  200  and  850 
mb;  4.  Temperature  difference  between  500  and  850  mb;  5.  Previous  24  hour  trend  of  700  mb 
equivalent  potential  temperature;  6.  200  mb  planetary  eddy  flux  convergence  at  700  km  radius;  7. 

Relative  angular  momentum  at  850  mb;  8.  Equivalent  potential  temperature  difference  between  700  and 
850  mb. 

Explained  Variance  of  Intensity  Change  (r2)  for  Petty  (1997) 


Table  2:  Summary  of  eastern  North  Pacific  tropical  cyclone  intensity  change 
forecasting  techniques  under  development  at  The  Ohio  State  University.  *  indicates  the 
five  forecast  variables  also  used  by  Petty  (1997). 
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CHAPTER  3 


DATA  AND  METHODOLOGY 

3.1  NHC  Postanalysis  Best  Track,  SST  and  ECMWF  Model  Data 

Best  track  data  indicating  postanalysis  determined  tropical  cyclone  positions  (0.1 
degree  increments)  and  intensities  (2.5  m  s'1  increments)  at  six-hour  intervals  during  each 
named  tropical  cyclone  (i.e.,  only  those  tropical  cyclones  that  reached  tropical  storm  or 
hurricane  intensity)  from  1992  to  1995  was  obtained  from  the  NHC 
(ftp://ftp.nhc.noaa.gov/pub/tracks/tracks.epa).  The  data  format  is  described  by  Davis  et 
al.  (1984).  The  positions  and  intensities  for  most  of  the  tropical  cyclones  in  the  best  track 
data  were  determined  by  satellite  techniques  (Dvorak  1990).  The  errors  associated  with 
using  the  Dvorak  intensity  specification  technique  are  not  specifically  known  for  the 
eastern  North  Pacific.  However,  in  a  study  of  North  Atlantic  tropical  cyclones.  Gaby  et 
al.  (1980)  found  the  mean  absolute  difference  (i.e.,  accuracy)  between  satellite-derived 
and  in  situ  intensity  measurements  was  4  m  s’1  and  the  mean  difference  (i.e.,  bias)  was  -2 
ms1.  Similarly,  in  a  study  of  western  North  Pacific  tropical  cyclones,  Martin  and  Gray 
(1993)  found  the  accuracy  of  satellite  estimates  was  10  m  s'1  and  the  bias  was  5  m  s’1. 

The  intensity  errors  in  the  eastern  North  Pacific  best  track  data  are  most  likely  on  the 
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order  of  the  errors  found  in  the  studies  cited  above. 

In  addition  to  the  best  track  data,  official  NHC  intensity  forecast  error  data  for 
these  tropical  cyclones  was  obtained  from  the  National  Climatic  Data  Center  (NCDC)  on 
3.5  inch  computer  diskettes  in  text  format. 

Climatology,  geography  and  persistence  variables  are  calculated  from  the  best 
track  data  in  the  same  manner  as  Hobgood  (1997).  The  only  exception  is  that  the 
calculation  of  potential  intensity  change  in  this  proposed  research  uses  real-time,  instead 
of  climatological,  SST  data.  Calculation  of  MPI  uses  SST  data  from  operational  analyses 
at  the  National  Centers  for  Environmental  Prediction  (NCEP).  This  SST  data  is  available 
from  a  data  archive  at  the  Columbia  University  in  binary  format 
(http://ingrid.ldgo.columbia.edu/descriptions/reynoldsweekly.html).  The  data  is 
optimally  interpolated  weekly  to  a  one-degree  grid.  Details  concerning  this  SST  data  are 
described  by  Reynolds  and  Smith  (1994).  SST  values  are  extracted  and  MPI  parameters 
calculated  using  the  same  computer  programs  from  Petty’s  (1997)  research. 

The  European  Centre  for  Medium-Range  Weather  Forecasts  (ECMWF)  NWP 
model  data  to  be  used  in  this  research  is  from  the  National  Center  for  Atmospheric 
Research  (NCAR)  data  archive  (http://www.scd.ucar.edu/dss/datasets/dsl  1  L2.html). 
Model  data  was  made  available  in  binary  format  to  the  Atmospheric  Sciences  Program  by 
the  Byrd  Polar  Research  Center  at  The  Ohio  State  University.  The  model  data  was 
optimally  interpolated  every  12  hours  (i.e.,  0000  and  1200  UTC)  to  a  2.5-degree  grid  with 
15  pressure  levels.  Synoptic-scale  environmental  forcing  predictors  identified  by  Petty 
(1997)  are  calculated  using  the  same  computer  programs  used  in  his  research.  The  use  of 
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ECMWF  model  data  to  evaluate  environmental  forcing  of  tropical  cyclones  was 
previously  shown  to  be  useful  by  Molinari  et  al.  (1992),  West  (1993a  and  b),  Rodgers  et 
al.  (1994a)  and  Rodgers  and  Pierce  (1995b). 

3.2  DMSP  Constellation  of  Satellites 

A  series  of  DMSP  satellites  (Block  5D-2)  have  been  launched  since  1987.  These 
satellites  include  flight  numbers  (F-#)  F-8,  F-9,  F-10,  F-l  1,  F-12  and  F-13.  All  of  these 
satellites  are  equipped  with  an  SSM/I  except  for  F-9.  The  SSM/I  on  F-12  failed  to 
function  properly  shortly  after  launch.  In  January,  1989,  the  F-8  experienced  a  partial 
failure  of  its  SSM/I  and  the  satellite  mission  was  terminated  in  August,  1991.  Satellites 
F-10  and  F-l  1  form  a  two  platform  constellation  for  observations  of  tropical  cyclones 
from  1992  to  1994.  The  additional  launch  of  F-13  in  1995  created  a  three  SSM/I 
constellation  for  the  1995  tropical  cyclone  season. 

Specific  information  about  the  orbits  of  the  relevant  DMSP  satellites  is  contained 
in  Table  3.  DMSP  satellite  orbits  are  near-circular,  near-polar  and  sun-synchronous  (i.e., 
the  satellite  crosses  the  Equator  at  the  same  local  time  each  day).  Inclination  of  the 
satellite  orbital  plane  is  98.8  degrees  from  the  Equator  resulting  in  a  retrograde  orbit.  The 
period  for  an  orbit  is  102  minutes  which  yields  14.1  orbits  per  day.  One  satellite,  F-10, 
did  not  reach  the  desired  orbit.  Therefore,  its  Equator  crossing  time  is  changing  at  a 
substantial  rate  shown  in  Table  3.  All  of  the  SSM/I-equipped  satellites  listed  in  Table  3 
are  currently  functioning  properly  and  within  specification  limits. 
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Flight- 

Number 

Launch  Date 

Equator 

Crossing 

Times  Upon 
Launch 
[Local  Standard 
Time] 

Equator 

Crossing 

Times  as  of  1 

Jan  97 

[Local  Standard 
Time] 

Maximum 

Altitude 

(km) 

Minimum 

Altitude 

(km) 

ascend 

descend 

ascend 

descend 

F-10 

1  Dec  90 

1942 

0742 

2223 

1023 

861 

726 

F-l  1 

28  Nov  91 

1702 

0502. 

1854 

0654 

878 

836 

F-13 

24  Mar  95 

1743 

0543 

1745 

0545 

877 

840 

Table  3:  DMSP  Orbital  Characteristics. 


3.3  SSM/I  Instrument  Specifications 

The  SSM/I  is  a  seven  channel,  four  frequency,  linearly  polarized,  passive 
microwave  radiometer  (Hollinger  et  al.  1990).  Spectral  channels  include  vertically  and 
horizontally  polarized  electromagnetic  radiation  at  19.4, 37.0  and  85.5  GHz.  Only 
vertically  polarized  radiation  is  measured  at  22.2  GHz.  Resolution  of  the  channels  is 
displayed  in  Table  4  (Hollinger  et  al.  1990;  Poe  1997;  Bauer  and  Grody  1995).  Absolute 
accuracy  of  the  SSM/I-measured  brightness  temperatures  is  less  than  3  K. 

Intercalibration  of  the  brightness  temperatures  measured  by  the  F-10,  F-l  1  and  F-13 
SSM/I’s  is  between  0.1  to  0.5  K  (Wentz  1995;  Poe  1997;  Colton  et  al.  1996). 
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Channel  Information 

Spectral 

Resolution 

Radiometric 

Resolution 

(K) 

Spatial  Resolution  (km) 

Sampling 

Resolution 

Center 

Frequency 

(GHz) 

Polarization 

(GHz) 

Along- 

Track 

Along- 

Scan 

(IFOV) 

Along- 

Scan 

(EFOV) 

(km) 

19.4±0.13 

Vertical 

0.24 

s  0.50 

69 

43 

70 

25.0 

19.4±0.13 

Horizontal 

0.24 

<;  0.48 

69 

43 

70 

25.0 

22.2±0.13 

Vertical 

0.24 

<;  0.55 

60 

40 

65 

25.0 

37.0±0.55 

|  Vertical 

0.90 

<;  0.37 

37 

28 

55 

25.0 

37.0±0.55 

Horizontal 

0.90 

s  0.37 

37 

— 

29 

55 

25.0 

85.5±0.80 

Vertical 

1.40 

s  0.58 

15 

13 

14 

12.5 

85.5±0.80 

Horizontal 

1.40 

<;  0.57 

15 

13 

14 

12.5 

Table  4:  SSM/I  Channel  Resolution  Information. 


SSM/I  scanning  characteristics  and  the  resulting  geometry  are  portrayed  in 
Figures  4  through  7.  The  SSM/I  is  a  conically  scanning  imager  with  an  active  scan  angle 
of  102.4  degrees  limited  by  the  body  of  the  spacecraft.  Depression  and  look  angles  of  45 
degrees  and  a  nominal  orbital  altitude  of  833  km  yield  an  incidence  or  zenith  angle  of 
53.1  degrees.  The  conical  scan,  active  scan  angle  and  orbital  altitude  result  in  a  ground 
swath  of  about  1394  km.  The  satellite  travels  with  a  speed  of  6.58  m  s'1  and  the  scan 
period  is  1.899  s.  Therefore,  the  satellite  ground  track  is  12.5  km  per  scan  which  nearly 
equals  the  spatial  resolution  of  the  85.5  GHz  channels.  At  a  sampling  interval  of  4.22  ms, 
128  samples  are  recorded  on  every  clockwise  (i.e.,  left-to-right)  scan  for  the  85.5  GHz 

channels.  The  remaining  channels  are  sampled  on  every  other  scan  at  a  8.44  ms  interval 
resulting  in  64  samples. 
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During  each  scan  the  SSM/I  is  calibrated  with  a  cold  target  (i.e.,  a  mirrored 
reflection  of  deep  space  which  has  a  brightness  temperature  of  about  2.7  K)  and  a  warm 
target  with  a  radiometric  temperature  near  300  K  obtained  by  three  precision  thermistors. 
The  conically  scanning  nature  of  the  SSM/I  yields  a  constant  incidence  angle  (i.e.,  53.1 
degrees)  and  sensor-to-ground  distance  throughout  each  scan.  This  constant  incidence 
angle  allows  for  invariable  planes  of  polarization  and  a  large  fixed  polarization 
difference.  Also,  the  constant  sensor-to-ground  distance  results  in  uniform  spatial 
resolutions  and  a  consistent  pathlength  along  each  scan.  Therefore,  the  limb  darkening 
problem  common  to  cross-track  scanners  is  avoided. 

3.4  SSM/I-Related  Data 

Data  measured  by  the  SSM/I  are  stored  onboard  the  satellite  until  transmitted  to 
ground  receiving  and  processing  stations.  During  initial  processing,  this  raw  data  is 
quality  controlled.  That  is,  the  data  are  decommutated,  deinterleaved,  bit-flipped  and 
restructured  into  properly  ordered  orbits.  Using  scan  calibration  data,  satellite  ephemeris 
and  attitude  corrections,  a  geolocated  antenna  temperature  file  called  a  Temperature  Data 
Record  (TDR)  is  created  containing  data  from  all  seven  channels.  The  geolocation  (i.e., 
ground  navigation  data)  accuracy  of  each  SSM/I  sample  is  less  than  6  km  (Hollinger 
1991).  A  correction  for  spill-over,  cross-polarization  coupling  and  antenna  pattern  effects 
caused  by  the  parabolic  reflector  used  in  the  SSM/I  is  required.  Therefore,  an  antenna 
pattern  correction  is  applied  to  the  TDR  to  produce  a  brightness  temperature,  Tu ,  file 
called  a  Sensor  Data  Record  (SDR).  Through  the  use  of  scientifically  developed 
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algorithms,  the  SDR  can  be  converted  into  a  file  containing  geophysical  parameter  data 
(e.g.,  rainfall)  called  an  Environmental  Data  Record  (EDR). 

SSM/I TDR  data  for  DMSP  satellite  flight  numbers  F-10  and  F-l  1  can  be 
obtained  for  the  1992  to  1995  tropical  cyclone  seasons.  In  addition,  F-l 3  data  are 
available  for  the  1995  season.  It  should  be  noted  that  due  to  SSM/I  problems  on  F-8  and 
F-l  2  during  the  period  of  interest,  data  from  these  specific  radiometers  is  not  available  for 
this  research.  The  specific  SSM/I  data  set  planned  for  use  in  this  research  is  maintained 
by  the  Naval  Research  Laboratory  in  Washington,  DC  and  Stennis  Space  Center, 
Mississippi.  Unfortunately,  it  is  not  a  contiguous  data  set.  That  is,  not  all  of  the  SSM/I 
orbits  from  1992  to  1995  are  contained  in  this  library. 

The  available  constellation  of  radiometers  virtually  assures  at  least  one 
coincidental  SSM/I  observation  during  a  tropical  cyclone  per  day  (Velden  et  al.  1989). 

The  probability  of  at  least  two  observations  per  day  is  greater  than  or  equal  to  0.89. 

Spatial  and  temporal  coverage  characteristics  of  SSM/I  orbits  are  graphically  summarized 
in  Figures  8  through  10. 

3.5  Query  Data  Library  for  SSM/I  Orbits  with  Coverage  of  Tropical  Cyclones 

Using  best  track  tropical  cyclone  positions  as  input,  a  computer  program 
developed  by  the  Naval  Research  Laboratory  in  Washington,  DC,  determined  orbits  with 
adequate  coincidental  coverage.  The  computer  program  used  radar-tracked  DMSP  orbital 
ephemeris  and  SSM/I  scanning  characteristics  to  perform  the  required  calculations. 

Output  of  this  program  was  input  to  a  Corel  Quattro  Pro  7  spreadsheet  database.  Using 
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database  functions,  orbits  were  selected  from  those  available  in  the  SSM/I  library  of  the 
Naval  Research  Laboratory  in  Washington,  DC  and  Stennis  Space  Center,  Mississippi. 

In  addition,  SSM/I  orbits  of  tropical  cyclones  located  west  of  140  degrees  West 
longitude  and  after  or  near  (i.e.,  within  55.5  km  of  a  coastline)  landfall  were  excluded 
from  consideration.  Determinations  of  distance  to  landfall  were  performed  using  the  best 
track  positions  and  coastline  positions  as  input  to  a  computer  program  used  at  the  NHC 
that  was  developed  by  Kaplan  (1992)  based  upon  previous  research  by  Merrill  (1987). 
Additional  cases  with  corresponding  best  track  intensities  less  than  15  ms*1  (i.e.,  30  knots 
or  Dvorak  (1990)  Current  Intensity  Number  2.0)  were  also  eliminated.  SSM/I-based 
tropical  cyclone  center  of  circulation  determinations  (see  Section  3.7  for  more  details)  are 
not  feasible  for  these  weaker  tropical  cyclone  intensities  (Hawkins,  personal 
communication,  1997;  Miller,  personal  communication,  1997). 

3.6  Process  SSM/I  Orbits  with  Coincidental  Coverage  of  Tropical  Cyclones 

A  computer  software  system  developed  by  the  Naval  Research  Laboratory  called 
TROPX  will  be  used  to  process  data  from  SSM/I  orbits  with  coincidental  coverage  of 
tropical  cyclones.  This  tropical  cyclone  satellite  data  processing  system  is  a  component 
of  the  Naval  Research  Laboratory  [NRL]  Satellite  Image  Processing  System  (NSIPS). 
TROPX  contains  display  and  analysis  software  tools  specifically  tailored  to  conduct 
multispectral  satellite  remote  sensing  studies  of  tropical  cyclones  (Helveston  et  al.  1996). 
TROPX  interfaces  with  Precision  Visuals-Workstation  Analysis  and  Visualization 
Environment  (PV-WAVE)  software  produced  by  Visual  Numerics  Incorporated.  The 
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combination  of  TROPX  and  PV~WAVE  allows  SSM/I  data  to  be  processed  from  the  raw 
data  format  found  in  the  Naval  Research  Laboratory^  archive  to  TDR,  SDR  and  EDR 
products. 

SSM/I  data  identified  as  containing  coincidental  coverage  will  be  extracted  from 
the  archive  and  transferred  to  a  disk  storage  device  attached  to  a  UNIX  workstation. 
Further  processing  and  database  functions  will  be  performed  by  TROPX  on  the  same 
workstation.  TROPX  can  be  controlled  locally  or  remotely  using  an  Internet  connection. 

3.7  Locate  the  Tropical  Cyclone  Center 

The  85.5  GHz  horizontally  polarized  channel  (i.e.,  an  SDR  product)  on  the  SSM/I 
provides  the  necessary  spatial  resolution  and  radiometric  response  to  accurately 
determine  the  center  of  tropical  cyclones  (Melton  1994;  Zhao  1994;  Alliss  et  al.  1993  and 
1992;  Rappaport  1991;  Glass  and  Felde  1989  and  1990;  Rhudy  1989;  Velden  et  al  1989). 
See  Figures  12  and  13  for  examples  of  85.5  GHz  horizontally  polarized  TB  images 
created  by  TROPX  that  are  useful  for  this  task.  The  latitude  and  longitude  of  the  center 
of  circulation  will  be  recorded  to  the  nearest  0.1  degree.  Along-track  scan  number,  along- 
scan  position  number  and  scan  time  of  the  center  of  circulation  will  also  be  recorded. 

This  information  and  interpolated  best  track  data  will  be  used  to  assess  the  accuracy  of 
tropical  cyclone  center  fixes  using  SSM/I  data. 
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3.8  Determine  if  SSM/I  Coverage  is  Adequate 

Once  an  accurate  tropical  cyclone  center  fix  is  accomplished,  only  an  orbit  with 
adequate  coverage  will  be  further  processed  from  a  brightness  temperature  (i.e.,  SDR)  to 
a  rainfall  (i.e.,  EDR)  image.  Swath  width  considerations  and  bad  or  missing  data  can 
prevent  adequate  coverage  by  a  coincidental  SSM/I  orbit.  Adequate  orbital  coverage  is 
defined  as  at  least  half  coverage  over  the  ocean  of  a  circular  area  within  444.4  km  of  a 
tropical  cyclone  center  111.1  km  or  more  from  land.  These  criteria  ensure  adequate 
rainfall  measurements  over  water  and  minimize  effects  due  to  landfall  (Merrill  1987). 

3.9  Apply  the  Operational  Rainfall  Algorithm 

The  rainfall  algorithm  (Ferraro  et  al.  1996)  used  in  this  study  was  selected  for  two 
reasons.  First,  it  has  the  capability  to  measure  both  convective  and  stratiform  rainfall 
over  the  ocean.  Second,  it  is  the  current  DoD  operational  rainfall  algorithm  (Colton  and 
Poe  1994).  See  Figures  14  and  15  for  examples  of  rain  rate  images  produced  by  the 
NESDIS/ORA  algorithm  using  TROPX. 

3.10  Examine  Rainfall  and  Intensity  Change  Correlations 

Using  Statistical  Package  for  Social  Sciences  (SPSS)  7.0  Base  Statistics  software, 
descriptive  statistics  and  data  analysis  tools  (e.g.,  scatter-diagrams,  histograms  and 
regression)  will  be  generated  (Glahn  1985).  The  average  and  time  tendency  of  average 
rainfall  parameters  for  regions  within  444.4  km  of  the  tropical  cyclone  center  will  be 
correlated  with  current  intensity,  previous  12-hour  intensity  trend  and  future  12-  to  72- 


41 


hour  intensity  changes  calculated  from  the  best  track  data.  Each  SSM/I  observation  will 
be  assigned  to  the  closest  best  track  time  interval  within  six  hours  of  0000,  0600,  1200  or 
1800  Uniform  Time  Coordinate  (UTC).  If  more  than  one  SSM/I  observation  corresponds 
to  the  same  record,  the  observation  closer  in  time  or  with  better  coverage  will  be  used. 

Circular  regions  are  designated  as  inner  core  (i.e.,  0  to  1 1 1.1  km),  core  (i.e.,  0  to 
222.2  km)  and  total  (i.e.,  0  to  444.4  km).  Annular  regions  are  the  outer  core  (i.e.,  1 1 1.1  to 
222.2  km)  and  exterior  (i.e.,  222.2  to  444.4  km).  These  regions  are  depicted  in  Figure  11. 
Average  rain  rates  for  these  regions  will  be  calculated.  Ratios  and  differences  of  the 
average  rain  rates  for  the  above  regions  will  also  be  calculated.  Additionally,  the  12-  and 
24-hour  time  tendency  of  the  rainfall  parameters  will  be  calculated.  The  resulting  number 
of  candidate  rainfall  parameters  is  5 1 .  These  parameters  are  summarized  in  Table  5 
below. 


Current 

12  hour  Tendency 
(Current  -  12  hours  prior) 

24  hour  Tendency 
(Current  -  24  hours  prior) 

Averages 

Ratios  of 
Averages 

Differences 
of  Averages 

Averages 

Ratios  of 
Averages 

Differences 
of  Averages 

Averages 

Ratios  of 
Averages 

Differences 
of  Averages 

inner 

outer 

core 

exterior 

total 

inner/outer 

inner/core 

inner/exterior 

inner/total 

outer/core 

outer/exterior 

outer/total 

core/exterior 

core/total 

exterior/totaJ 

inner-outer 

core-exterior 

inner 

outer 

core 

exterior 

total 

inner/outer 

inner/core 

inner/exterior 

inner/total 

outer/core 

outer/exterior 

outer/total 

core/exterior 

core/total 

exterior/total 

inner-outer 

core-exterior 

inner 

outer 

core 

exterior 

total 

inner/outer 

inner/core 

inner/exterior 

inner/total 

outer/core 

outer/exterior 

outer/total 

core/exterior 

core/total 

exterior/total 

inner-outer 

core-exterior 

Table  5:  Candidate  Rainfall  Parameters. 
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Stratifications  will  be  used  to  determine  possible  factors  affecting  significant 
rainfall  parameter  and  intensity  change  correlations.  These  stratifications  include:  current 
intensity,  rate  of  intensity  change,  month,  location,  SST,  vertical  shear  of  the  horizontal 
wind,  translation  speed  and  landfall.  Intensity  stratifications  are  tropical  depression  (i.e., 
less  than  17  m  s'1),  tropical  storm  (i.e.,  17  to  32  m  s'1),  hurricane  (i.e.,  33  to  50  m  s'1)  or 
strong  hurricane  (i.e.,  greater  than  50  m  s'1).  Rate  of  intensity  change  classes  are  12-hour 
change  in  intensity  greater  than  1 1  m  s'1  (i.e.,  rapid  strengthening),  1 1  to  -1 1  m  s*1  (i.e., 
normal)  or  less  than  -11ms'1  (i.e.,  rapid  weakening).  Month  stratifications  are  peak  of 
season  (i.e.,  August  and  September),  early  season  (i.e.,  before  August)  and  late  season 
(i.e.,  after  September).  Location  stratifications  are  north  or  south  of  18  degrees  North 
latitude  and  east  or  west  of  1 10  degrees  West  longitude.  SST  stratifications  are  less  than 
26  °C  or  greater  than  or  equal  to  26  °C.  Vertical  shear  of  the  zonal  component  of  the 
horizontal  wind  between  200  and  850  mb  stratifications  are  weak  (i.e.,  less  than  15  m  s'1) 
or  strong  (i.e.,  greater  than  or  equal  to  15  m  s'1).  The  stratifications  for  12-hour  average 
translation  speed  are  slow  (i.e.,  less  than  or  equal  to  10  m  s'1)  or  fast  (i.e.,  greater  than  10 
m  s'1).  Landfall  stratifications  are  whether  or  not  the  tropical  cyclone  made  landfall 
during  its  existence. 

3.11  Develop  Rainfall  Parameter  and  Intensity  Change  Prediction  Model 

Those  candidate  rainfall  parameters  found  to  be  significantly  correlated  with 
intensity  changes  at  the  0.95  level  will  be  rank  ordered  by  correlation  as  potential 
predictors  for  multiple  linear  regression  model  development.  Rainfall  parameters  highly 
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correlated  with  12-,  24-,  36-,  48-,  60-  and  72-hour  intensity  changes  will  be  selected  for 
future  analysis.  Also,  selected  parameters  will  be  examined  to  ensure  that  they  are  not 
highly  correlated  with  each  other.  The  selected  potential  predictors  (i.e.,  independent 
variables)  and  dependent  variables  (i.e.,  12-  to  72-hour  intensity  changes)  for  the  model 
development  data  set  (i.e.,  1992  to  1994)  will  be  normalized.  That  is,  all  variables  will 
have  the  development  data  set  mean  for  that  variable  subtracted  from  it  and  this 
difference  divided  by  the  standard  deviation  for  the  development  data  set.  This  allows  for 
dimensionless  comparison  of  the  resulting  regression  coefficients  between  predictors  and 
across  dependent  variables  (i.e.,  intensity  changes  for  the  different  forecast  time  periods). 
Via  the  SPSS  step-wise  multiple  linear  regression  model  development  tool,  a  prediction 
model  based  upon  rainfall  parameters  will  be  constructed.  The  statistical  significance 
level  for  entry  into  the  regression  model  will  be  0.95  and  0.90  to  remain. 

3.12  Combine  Rainfall  Parameters  with  Hobgood  and  Petty  Prediction  Models 

To  evaluate  the  impact  of  rainfall  parameters  on  other  intensity  change  prediction 
models  under  development  for  the  eastern  North  Pacific  (Hobgood  1997;  Petty  1997), 
combined  multiple  linear  regression  models  will  be  constructed  in  a  similar  manner  as  the 
rainfall  parameter  model.  Predictors  described  by  Hobgood  (1997)  will  be  computed  for 
the  six-hour  interval  best  track  data  records  with  corresponding  rainfall  parameters 
selected  using  the  same  model  development  data  set  as  Section  3.11.  Hobgood’s  (1997) 
climatology ,  geography  and  persistence  predictors  and  rainfall  parameter  predictors 
selected  in  Section  3.1 1  will  be  input  as  predictor  candidates  into  the  SPSS  multiple 
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linear  regression  model  development  tool  in  the  same  manner  as  described  in  Section 
3.1 1.  To  combine  the  rainfall  parameter  predictors  with  Petty’s  (1997)  model,  the  rainfall 
parameters  must  correspond  to  0000  or  1200  UTC.  Those  rainfall  parameters  derived 
from  F-l  1  (and  F-13)  SSM/I  orbits  will  normally  be  within  three  hours  of  either  0000  or 
1200  UTC.  However,  F-10  orbits  will  generally  be  closest  to  0600  or  1800  UTC.  To 
prevent  large  data  gaps,  F-10  rainfall  parameter  observations  can  be  used  for  the  nearest 
time  interval  within  six  hours  of  0000  or  1200  UTC,  if  a  corresponding  F-l  1  observation 
is  missing.  Petty  s  (1997)  climatology,  persistence  and  synoptic-scale  environmental 
parameter  predictors  and  rainfall  parameter  predictors  will  be  used  to  develop  a  third 
multiple  linear  regression  model  as  described  above. 

3.13  Evaluate  Intensity  Change  Prediction  Models 

The  skill  of  the  three  models  developed  in  Sections  3.1 1  and  3.12  will  be 
evaluated  by  examining  the  intensity  change  forecast  errors  they  produce  on  a  data  set 
(i.e.,  1995)  that  is  independent  of  the  data  set  used  to  create  the  models.  Also,  the  errors 
will  be  compared  to  those  produced  by  the  following  methods:  Hobgood  (1997),  Petty 
(1997),  official  NHC  forecasts,  simple  persistence  of  the  past  12-hour  intensity  change 
trend  and  predicting  no  change.  Two  types  of  forecast  errors  will  be  examined  for  each 
forecast  interval  from  12  to  72  hours.  Mean  intensity  change  forecast  error  (i.e.,  bias  in 
units  of  m  s  ')  is  the  average  value  of  the  difference  between  the  forecasted  intensity 
change  and  the  intensity  change  that  actually  occurred.  Mean  absolute  intensity  change 
forecast  error  (i.e.,  accuracy  in  units  of  m  s'1)  is  the  average  absolute  value  of  the 
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difference  between  the  forecasted  intensity  change  and  the  intensity  change  that  actually 
occurred. 

For  the  regression-based  models,  the  predictor  variables  for  the  evaluation  data  set 
(i.e.,  1995)  will  be  normalized  using  the  mean  and  standard  deviation  from  the  model 
development  data  set  (i.e.,  1992  to  1994).  These  normalized  values  will  be  input  to  the 
regression  models.  The  output  will  be  converted  from  normalized  values  to  intensity 
changes  for  comparison  purposes.  The  forecast  errors  described  above  will  be  calculated 
for  all  the  methods.  Both  types  of  forecast  errors  for  each  method  will  be  graphically 
displayed  using  Quattro  Pro.  To  understand  how  well  each  method  performs  in  various 
situations,  forecast  errors  will  also  be  stratified  by  current  intensity  and  rate  of  intensity 
change  categories  described  in  Section  3.10.  The  statistical  significance  of  differences  in 
errors  produced  by  the  different  methods  will  be  assessed  by  performing  paired  t  tests  at 
the  0.95  level. 

3.14  Conduct  Case  Studies 

Tropical  cyclones  with  particularly  good  spatial  and  temporal  coverage  by  SSM/I 
orbits  will  be  candidates  for  in-depth  case  studies.  Also,  strong  hurricanes  with 
maximum  intensities  greater  than  50  m  s'1  are  significant  cases  that  might  prove  to  be 
insightful  upon  close  examination.  Additionally,  those  cases  exhibiting  rapid  intensity 
changes  or  the  concentric  eyewall  cycle  could  provide  interesting  examples  of  forecast 
model  skill  and  operational  implementation  issues  for  challenging  situations.  Tropical 
cyclones  Tina  (1992)  and  Olivia  (1994)  were  subjects  of  aerial  reconnaissance  research 
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missions  and  would  be  excellent  candidates  for  case  studies. 


3.15  Survey  Rainfall  Characteristics  of  Eastern  North  Pacific  Tropical  Cyclones 

Averages  and  standard  deviations  of  rain  rates  for  the  areas  shown  in  Figure  1 1 
will  be  computed  for  the  intensity  and  intensity  change  categories  in  Section  3.10. 
Calculations  will  use  SSM/I  observations  from  1992  to  1995.  These  calculations  will  be 
further  stratified  by  time  of  day  in  an  effort  to  detect  a  diurnal  cycle  of  rainfall. 

3.16  Analyze  Tropical  Cyclone  Center  Position  Errors  Using  the  SSM/I 

The  average  and  standard  deviation  of  the  distance  errors  between  SSM/I 
determined  tropical  cyclone  centers  of  circulation  (see  Section  3.7)  and  interpolated  best 
track  locations  will  be  calculated  for  1992  to  1995.  Errors  will  be  stratified  by  current 
intensity  categories  in  Section  3.10. 
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CHAPTER  4 


ACCOMPLISHMENTS  AND  PRELIMINARY  RESULTS 


Time  Period 

Number  of 
tropical  cyclones 
with  maximum 
intensity  of 
tropical  storm 

Number  of 
tropical  cyclones 
with  maximum 
intensity  of 
hurricane 

Total  number  of 
tropical 

cyclones  for  the 
time  period 

Number  of 
tropical  cyclones 
with  at  least  one 
coincidental 

SSM/I  orbit 

Number  of 
coincidental 

SSM/I  orbits 
identified  for 
further  processing 

1992“ 

10 

14 

24 

23 

207 

1993b 

4 

10 

14 

12 

93 

1994“ 

8 

9 

17 

6 

49 

1992-1994 

22 

33 

55 

41 

349 

1995d 

3 

7 

10 

10 

127 

1992-1995 

25 

40 

65 

51 

476 

Climatology 

(1966-1995)“ 

7 

9 

16 

l 

l 

l 

l 

1 

Table  6:  SSM/I  coverage  of  eastern  North  Pacific  tropical  cyclones  and  number  of  orbits 
selected  for  further  processing.  Climatology  values  provided  for  comparison  (“Lawrence 
and  Rappaport  1994;  bAvila  and  Mayfield  1995;  “Pasch  and  Mayfield  1996;  dAvila  and 
Rappaport  1996;  “Rappaport  and  Mayfield  1997). 


A  database  of  SSM/I  orbits  with  coincidental  coverage  of  eastern  North  Pacific 


tropical  cyclones  has  been  created  for  the  seasons  of  1992  to  1995.  476  orbits  have  been 
identified  for  further  processing  at  the  Naval  Research  Laboratory  Marine  Meteorology 
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Division  in  Monterey,  California.  The  dependent  data  set  (i.e.,1992  to  1994)  to  be  used 
in  model  development  contains  349  orbits.  There  are  127  orbits  in  the  independent  data 
set  (i.e.,  1995)  to  be  used  for  model  evaluation.  Information  about  coincidental  coverage 
of  SSM/I  orbits  selected  for  further  processing  is  summarized  in  Table  6  above. 

Some  of  the  climatology,  geography,  persistence  and  synoptic-scale 
environmental  variables  required  by  the  intensity  change  prediction  models  of  Hobgood 
(1997)  and  Petty  (1997)  have  been  calculated  for  1993  to  1995.  The  main  focus  of 
current  research  activities  is  extraction  of  SSM/I  orbits  from  the  archive,  processing  of 
the  SSM/I  data,  generation  of  the  85.5  GHz  horizontally  polarized  TB  images  for  center 
fixing,  application  of  the  NESDIS/ORA  rainfall  algorithm  and  calculation  of  rainfall 
parameters.  Tasks  related  to  the  SSM/I  data  will  require  the  author  to  travel  to  the  Naval 
Research  Laboratory  in  Monterey,  California. 

It  is  expected  that  rainfall  parameters  will  show  predictive  skill  in  determining 
intensity  change  and  aid  efforts  to  develop  a  better  objective  guidance  product  for  the 
intensity  prediction  problem.  A  summary  of  this  proposed  research  was  recently 
presented  by  the  author  at  the  22nd  Conference  on  Hurricanes  and  Tropical  Meteorology 
sponsored  by  the  American  Meteorological  Society  (West  1997).  The  author  expects  to 
defend  his  dissertation  by  the  Summer  of  1998. 
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INFRARED  COOLING 


Figure  1:  Interactions  of  synoptic-scale  water  and  motion  fields  (Atlas  and  Thiele  1981). 
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Figure  2:  Interactions  involving  microwave  radiation  (Rao  et  al.  1990). 
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Surface  £~0.5  for  water 
£-0.9  for  land 


Figure  4:  Conical  scanning  radiometer  (Skou  1989). 
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Figure  5:  SSM/I  scan  geometry  (Hollinger  et  al.  1990). 
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Figure  6:  SSM/I  scan  characteristics  (Spencer  et  al.  1989). 


Along  Track 
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Figure  7:  SSM/I  spatial  sampling  illustrating  along-scan  changes  (Hollinger  1 989). 


•e  8:  Equatorial  view  of  successive  SSM/I  orbits  with  swath  widths  (Hollinger  1989). 


Figure  9:  Polar  view  of  successive  SSM/I  orbits  with  swath  widths  (Hollinger  1989). 
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1400  km  Swath 


Figure  10:  Coverage  by  one  SSM/I  in  24  hours.  Dark  areas 
(Hollinger  1989). 
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Total 

Figure  11:  Circular  and  annular  regions  to  be  used  for  average  rain  rate  calculations. 
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Figure  12:  15-panel  85.5  GHz  horizontally  polarized  Ts  images  created  by  TROPX.  Images  depict  Tropical  Cyclone  Tina 
(1992)  from  20-28  September.  Number  in  the  upper-left  comer  of  each  image  is  interpolated  best  track  intensity  in  knots. 
Date,  time  and  DMSP  flight  number  are  indicated  below  each  image. 
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Figure  13:  15-panel  85.5  GHz  horizontally  polarized  TB  images  created  by  TROPX.  Images  depict  Tropical  Cyclone 
Tina  (1992)  from  28  September  to  10  October.  Number  in  the  upper-left  comer  of  each  image  is  interpolated  best  track 
intensity  in  knots.  Date,  time  and  DMSP  flight  number  are  indicated  below  each  image. 
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Figure  14:  15-panel  NESDIS/ORA  rainfall  algorithm  images  created  by  TROPX.  Images  depict  Tropical  Cyclone  Tina 
(1992)  from  20-28  September.  Number  in  the  upper-left  comer  of  each  image  is  interpolated  best  track  intensity  in 
knots.  Date,  time  and  DMSP  flight  number  are  indicated  below  each  image. 
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Figure  15:  15-panel  NESDIS/ORA  rainfall  algorithm  images  created  by  TROPX.  Images  depict  Tropical  Cyclone  Tina 
(1992)  from  28  September  to  10  October.  Number  in  the  upper-left  comer  of  each  image  is  interpolated  best  track 
intensity  in  knots.  Date,  time  and  DMSP  flight  number  are  indicated  below  each  image. 


